Dark energy exponential potential models as curvature quintessence 
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It has been recently shown that, under some general conditions, it is always possible to find a 
fourth order gravity theory capable of reproducing the same dynamics of a given dark energy model. 
Here, we discuss this approach for a dark energy model with a scalar field evolving under the action 
of an exponential potential. In absence of matter, such a potential can be recovered from a fourth 
order theory via a conformal transformation. Including the matter term, the function f[R) entering 
the generalized gravity Lagrangian can be reconstructed according to the dark energy model. 
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I. INTRODUCTION 

The era of precision cosmology has opened the way to 
a new picture of the universe depicted as a spatially flat 
manifold with a subcritical matter content undergoing 
a phase of accelerated expansion. Although totally un- 
expected only few years ago, this scenario is nowadays 
fully accepted as the standard one, thanks to the im- 
pressive amount of astrophysical evidences on different 
scales. The anisotropy spectrum of cosmic microwave 
background radiation jj, , the SNela Hubble diagram 
HI, the large scale structure and the matter power 
spectrum determined by the Lya forest data ^ repre- 
sent observational cornerstones that put on firm grounds 
the picture of this kind of universe model. 

While observationally well founded, this scenario is 
plagued by the unsolved puzzle of dark energy, the neg- 
ative pressure fluid which seems to be necessary to both 
close the universe and drive its accelerated expansion. 
Investigating the nature and the properties of this mys- 
terious component is one of the most exciting problems of 
modern cosmology, with proposals ranging from the clas- 
sical cosmological constant Q , to scalar field quintessence 
and unified dark energy models in which a single fluid 
with an exotic equation of state accounts for both dark 
energy and dark matter 

Notwithstanding the strong efforts made to solve this 
puzzle, none of the proposed explanations is fully satis- 
factory and free of problems. This disturbing situation 
has motivated much interest towards radically different 
approaches to the problem of dark energy. After all, dark 
energy is usually introduced only as a tool to explain the 
cosmic acceleration, while there is no direct evidence of 
its existence as a separate fluid. Motivated by this con- 
sideration, it has been suggested that cosmic speed up 
could be an evidence for the need of new physics rather 
than a new fluid. Much interest has then been devoted. 



for example, to models according to which standard mat- 
ter is the only physical ingredient, while the Friedmann 
equations are modified as in the Cardassian exp ansion 
, the Dvali - Gabadadze - Porrati scenario J13| or the 
improved renormalization group approach . 

In this same framework, higher order theories of grav- 
ity also represent a valid alternative to the dark energy 
approach. In these models, also referred to as curvature 
quintessence |l4l |. the gravity Lagrangian is generalized 
by replacing the Ricci scalar curvature R with a generic 
function f{R), so that an effective dark energy- like fluid 
appears in the Friedmann equations and drives the ac- 
celerated expansion. Reasons for considering non linear 
curvature terms in the gravity Lagrangian can be found, 
for instance, in quantum effects on curved spacetimes or 
in certain low energy limits of string/M theories. More- 
over, early time inflationary behaviours can be naturally 
recovered from nonlinear gravity Lagrangians thus fur- 
ther motivating to study possible effects of such terms 
on the late time universe expansion. Different models of 
this kind have been explored and tested against observa- 
tions also considering the two possible formulations that 
are obtained adopting the metric 0, 0, 0, or the 
Palatini [ll[ll|23| formulations. 

Recently, it has been shown ^1] that, under some quite 
general conditions, the dynamics (i.e., the scale factor 
and the expansion rate) of a given dark energy model 
may be, in principle, reproduced as a byproduct of a f(R) 
theory of gravity^. The method developed in j23| makes 
it possible to derive the f{R) theory which reproduces a 
given H{z) and represents therefore a sort of bridge be- 
tween two physically different scenarios. Here, we apply 
the method to a quintessence model where the evolution 
of the scalar fleld is determined by an exponential poten- 
tial [2^ 13 . Since this model has been shown to fit well 
a large set of observational constraints [2^ |23| , we 
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^ See also |22| for some comments on the uniqueness of the recon- 
struction in the matter only case. 
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are thus able to recover a f{R) theory that reproduces 
the same data set with the same accuracy, so being mo- 
tivated a posteriori and showing that the two approaches 
are phenomenologicaUy equivalent. 

The plan of the paper is as follows. In Sects. II and III, 
we shortly review the main features of the dark energy 
model and of the curvature quintessence scenario, respec- 
tively. In Sect. Ill, the procedure adopted to recover the 
f{R) theory representing the higher order counterpart of 
a given dark energy model is presented. In Sect. IV, we 
use conformal transformations to show the equivalence 
between a class of f{R) theories and a given scalar field 
model with exponential potential. Since this equivalence 
only holds in absence of matter. Sect. IV presents the re- 
sults obtained applying the method outlined in Sect. Ill 
to recover the f{R) Lagrangian corresponding to the ex- 
ponential potential model in presence of the matter term. 
Some caveats related to our technique are commented 
upon in Sect. V, where we also summarize and conclude. 



II. THE EXPONENTIAL POTENTIAL DARK 
ENERGY MODEL 

A standard recipe to explain the observed cosmic accel- 
eration consists in considering a scalar field ip minimally 
coupled to gravity yielding a homogenously distributed 
fluid with a negative equation of state. The energy den- 
sity and pressure for such a fluid read^ 



H = 



2 [2 {3Ho - 2) r2 + 4 - 3Ho\ 
3r [2 (3Ho - 2) t2 + 4 - 3Hq] 



(4) 



with T — t/to, to the present age of the universe and Tio 
the Hubble constant in units of to and hereon the sub- 
script will refer to quantities evaluated today. Taking 
to as unit of time, the model is fully characterized by 
the parameter Tio- In particular, the present day matter 
density parameter turns out to be 
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(5) 



which, for the best fit value Tio = 0.97 |27|, gives VIm 
0.26, in very good agreement with recent estimates. 



III. CURVATURE QUINTESSENCE 

Much interest has been recently devoted to the so 
called curvature quintessence according to which the uni- 
verse is filled by pressureless dust matter only and the ac- 
celeration is the result of the modified Friedmann equa- 
tions obtained by replacing the Ricci scalar curvature R 
with a generic function /(i?) in the gravity Lag rangian. 
The Friedmann equations therefore read [Tj, [I3 : 



Pv 



V{ip) 



(1) 



with an overdot denoting derivation with respect to cos- 
mic time t and V{ip) the self interaction potential. As it 
is clear, the choice of V{ip) plays a key role in determining 
the dynamics of the universe. 

As an interesting a ppli cation, we consider the simple 
exponential potential |23L I2II25I I2ll27j : 



= -w. 



curvHcurv i 



(6) 



(7) 



where the prime denotes derivative with respect to R, 
Pcurv is the energy density of an effective curvature fluid 
given as : 



V{ip) oc exp 




(2) 



which presents the nice feature that the corresponding 
Friedmann equations may be analytically solved. We re- 
fer the interested reader to [23| , while here we only report 
the quantities we will need in the following. In a spatially 
flat universe, the scale factor a and the Hubble parameter 
H = a/a read (HH^ : 



= [{SHo - 2y + 4 - 3Ho] 



(3) 



1 fl 



f'{R) 12 



[f{R) - Rf'iR)] - 3HRf"{R) , 



and the barotropic factor of the curvature fluid is : 



= -1 



Rf'iR) + R R!"'{R) - Hf"{R) 



[f{R) ~ Rf'iR)] /2 - iHRf'iR) 



(8) 



(9) 



Assuming that there is no interaction between matter 
and curvature terms, the continuity equation for pcurv 
reads El : 



We use here natural units such that SttGjs 



'iHi\ + Wcurv)p c 



SH^nMRfiR) _ 



[riR)r 



a-' (10) 
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which is identically satisfied as can be easily shown using 
Eq. (|SJ) and expressing the scalar curvature R as function 
of the Hubble parameter : 



R^-6{H + 2H 



(11) 



Combining Eqs.® with Eq.© and using the definition 
of H, one finally gets the following master equation for 
the Hubble parameter [23| : 



H = - 



2f\R) 



[sH^^Ma-^ + Rf"{R)- 



+R 



RriR)-Hr{R)]} 



(12) 



Inserting Ea. l|ll|l into Ea. H12() . one ends with a fourth 
order nonlinear differential equation for the scale factor 
a{t) that cannot be analitically solved also for the sim- 
plest cases (for instance, f{R) oc i?" unless dust matter 
contribution is discarded). Moreover, although techni- 
cally feasible, a numerical solution of Ea. H12|) is plagued 
by the large uncertainties on the boundary conditions 
(i.e., the present day values of the scale factor and its 
derivatives up to the third order) that have to be set to 
find out a{t) by solving Ea.fT^. 

Given these mathematical difficulties, a different ap- 
proach has been proposed in (hereafter CCT) where 
Ea. (|12|l is considered as a way to determine f{R) rather 
than a{t). Rearranging the different terms suitably, CCT 
obtained a linear third order differential equation for / in 
terms of the redshift z = 1/a — 1 (having set qq = 1) that 
can be easily solved numerically for a given H{z). By 
this method, it is then possible to find out which f{R) 
theory reproduces the same dynamics of a given dark en- 
ergy model, thus showing a formal equivalence between 
these two radically different approaches. 

CCT developed the method using the redshift z as in- 
tegration variable since it is common to have an analyti- 
cal expression for the Hubble parameter as function of z 
rather than t. However, this is not the case for the mod- 
els considered here, so that we rewrite the main formulae 
in CCT using t as integration variable. The equation 
determining f{t) is then'^ : 



.11 

dt^ 



^3(i)^+W2(t)^+7^l(t) 



dt^ 



dt 



with : 



Hi = 2HR + HR + 2R^R-^ - R 



-3H^nMR^a-^t) , 
(13) 

(14) 



2R + HR 



(15) 



With an abuse of notation, we write f(t) rather than f[R(t)]. 



~ R ■ 



(16) 



where R is given by Ea. Hll() . In order to integrate 
Eg. lfT^ . we need to specify boundary conditions that 
are more conveniently assigned at the present time. We 
slightly generalize here the discussion presented in CCT. 
First, let us remember that, in a fourth order theory, we 
may define an effective gravitational constant as Gg/ / = 
Gn / f'{R), with Gn the usual Newtonian gravitational 
constant. Its rate of variation will be given as : 



G 



G 



eff 



1 f"{R) dR 
h f'{R) ^ 



(17) 



It is quite natural to assume that the effective and the 
Newtonian gravitational couplings take the same values 
today, so that we get the condition : 



1 



f'{Ro) 

Evaluating Ea. l(T7| for t — to (i.e., r 



(18) 



1), we may deter- 
mine /"(i?o) provided an estimate of {Gef f /Gef f)t=tQ 
is given. Since in our theory Gn is constant, we may 
assume that the measurements of the variation of Gn 
[23 actually refers to Gejf and use these results to get 
an estimate of {Gef f /Gef f)t=to- We thus take as our 
boundary condition : 



/"(i?o) 



G 



eff 



G 



eff 



obs 



dR 
d7 



(19) 



having used Ea. lfH^ . Finally, inserting Eas. ltT^ and itTi^ 
into Eq. (|SJ) and then in © evaluated today, we get : 



fiRo) 



6H^ i^l - Qm 
'dR^ 



QHo 



dt 



Rq 
/"(i?o) 



f{Ro 



(20) 



It is then straightforward to show that the following 
boundary conditions descend from Eas. lfT?^ - : 



dt^ 



t = to 



t = to 

dR 
'dt 



dR 
'dt 



/'(i?o) 



(21) 



f"{Ro) + 



t=to 



fR 
dt^ 



t=to 



fiRo) 

(22) 

that have to be used, together with Ea. H20|) . to numer- 
ically solve Ea. H13|) . Combining the solution thus ob- 
tained for f{t) with R{t) evaluated through Ea. (|lll) . one 
finally finds f{R) thus recovering the higher order theory 
that mimicks the assigned dark energy model. We refer 
the reader to [2lj | for some interesting examples. 
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IV. CONFORMAL EQUIVALENCE 

The procedure described above may be easily imple- 
mented in the case of the exponential potential model 
presented in Sect. II. Neglecting the matter term (i.e., 
setting = 0), it is however possible to show that 
there exists a more intimate relation between /(i?) the- 
ories and a scalar field evolving under the action of an 
exponential potential. 

To this aim, let us remember that, in general, given a 
/(i?) theory in the Jordan frame, the conformal trans- 
formation 

— o2 



with Vl = exp {lo) — \J f'{R), leads, in the Einstein frame, 
to a minimally coupled scalar field theory whose self in- 
teraction potential is given by : 



having used df/d(p = {df /dR){dR/d(p) = f{R)dR/dip. 
Now, it can be easily shown that posing 



dV 

dip 



-XV , 



(28) 



with A an arbitrary constant, leads to / = /o^"i where 
/o is an integration constant and n = (AK — X) / {2K — X) . 
Therefore, starting now from 



f{R) = i?" , 

the potential in Ea. (|25|l becomes 

1 



(29) 



V = V{R) 



2n2 



fjR) - Rf'jR) 



(23) 



where tilted quantities refer to the Einstein frame. Let 
us now consider the particular choice : 



1 — rt 

V = V{lp) = —-\n\ exp(-A¥5) , 



with A = 2K{n - 2)/{n - 1). 
Assuming K = -y/l/G and A 



(30) 



fiR)^R+- 



(24) 



and recall that the conformal factor may also be rescaled 
as a; = Kip with K an (up to now) arbitrary constant 
[30l |. Imposing f'{R) = cxp (w) and solving with respect 
to R, the potential H23|l may be rewritten as : 



V{cp) 



■ exp 



(25) 



where we have assumed a < and hereon we will use 
nontilted quantities for simplicity. Choosing K = 1/6 
and considering the limit exp {2K(p) >> 1, the potential 
()25|l asymptotically behaves as : 



V{(p) - exp I -xl-ifi 



which is just the exponential potential model we are con- 
sidering in Eq.Q. 

On the other hand, starting directly with f'{R) = 
exp {2Klp), the potential in Ea. l|25(l is soon written as 



V = 



f -Rcxp {2Klp) 
2 exp {-iKtp) 



(26) 



so that 



^ ^ ^ RfiR) - 2f{R) 



dip 



(27) 



-1 



(31) 



so that we get the nice result that a higher order theory 
of gravity (without matter, anyway) with 



f{R) 



R 



(32) 



is conformally transformed into a usual minimally cou- 
pled general relativistic theory with a massive scalar field 
p only, whose potential is 



V{p:) = exp 




(33) 



By means of conformal transformations, it is then possi- 
ble to find out an intimate link between two different sce- 
narios, namely a particularly simple class of /(i?) theories 
and a minimally coupled scalar field with specific expo- 
nential potential. It is worth stressing, however, that the 
above results only hold in the absence of matter. While 
it is still possible to work out a conformal transformation 
also in presence of matter, the self interaction potential 
may be no more expressed as in Eg . 1)23(1 . Although such 
a problem may be in principle escaped by resorting to nu- 
merical methods, it is nevertheless problematic to assess 
the meaning of conformal transformations, since there is 
a still open debate regarding the physical equivalence of 
two conformally (and hence mathematically) equivalent 
theories (see [3fl] and references therein). 
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FIG. 1: Reconstructed f{t) as function of r = t/to. For 
clarity, we plot If = log(///o) rather than / itself. The 
upper left panel refers to the period near the recombination 
epoch, while the lower right panel extends in the near future. 



V. f(R) AND EXPONENTIAL POTENTIAL 

The dark energy model presented in Sect. II has nice 
features both from the mathematical and the physical 
points of view. Indeed, the main dynamical quantities 
(such as the scale factor and the expansion rate) may 
be written analytically, while it very well matches the 
SNela Hubble diagram, the growth index and the posi- 
tion of the peaks in the CMBR anisotropy spectrum. It 
is, therefore, interesting to apply the method outlined 
above to find out the f{R) counterpart of this model. 
To this aim, we just need to insert the expressions for 
a{t) and H{t) into Eas. H13l) - I21|) and numerically deter- 
mine first /(<) and then f{R). As a preliminary step, 
we set 7^0 = 0.97 as obtained in f27!|. To estimate 
/"(.Ro) through Eq.lSI), we use to = 13.7 Gyr 2] and 
{Geff /Geff)obs ~ 1.6x10^^^ yr~^ that is the upper limit 
determined through heliosismological studies 3lJ. Note 
that this choice is somewhat arbitrary since there are 
a lot of different estimates obtained with radically dif- 
ferent methods For instance, using Lunar Laser 
Ranging, Williams et al. ^| have recently obtained 
Geff /Geff = (4 ± 9) X 10^3 yj.-i. Fortunately, most 
of the estimates available in literature are all consistent 
within the uncertainties so that we are confident that any 
bias is not introduced in the reconstruction process. 

Before discussing the results, there is a conceptual 
point to clarify. The method described in Sect. Ill rep- 
resents a sort of bridge between two different scenarios. 
In particular, their matter contents could be different, so 
that we should define both fi^"''^ and fl^Jj to denote this 
quantity in the two different models. In principle, there is 
no reason why il^"''" — H.'^/f should hold. However, since 

fi^{ is close to the fiducial value (f^M ~ 0.3) suggested by 
model independent estimates (e.g., from galaxy clusters 
abundance), we take flf}"^^ = ri^{. 

The reconstructed f{t) for the exponential potential 
model is shown in Fig.^ where we plot log(///o) as 



FIG. 2: Reconstructed f{R) for the exponential potential 
model. In both panels, rj = {f / fo)/{R/ Ro) is plotted for 
sake of clarity. Left panel shows the result of the reconstruc- 
tion process over the full time range 10~^ < r < 1.5, so that 
we report log Rs (with Ra = R/Rq) on the abscissa for sake of 
clarity. In the right panel, we zoom on the epoch 0.1 < r < 1 
corresponding to < « < 4.7. Note that 7? is a decreasing 
function of t so that the time increases towards right. 



function of the dimensionless time r. It is, however, more 
interesting to look at Fig.|21 where r] = (///o)/(i?/i?o) 
is plotted in terms of the scalar curvature R normalized 
to its present day value Rq. Although an analytical so- 
lution has not been obtained, we have found that the 
numerically reconstructed Lagrangian may be very well 
approximated^ as : 



f{R) ^ vsfo 



X O 



R 

Vi 1 



(34) 



R 




where the parameters (77^, 771, 772, a, /3) depend on the 
value adopted for /"(i?o). For our case, it is: 

{ri,,rii,r]2,a,(3) = (1.07,0.029,0.810,0.232,0.221) . 

The approximating function in Eq.(33) indeed shows to 
be a useful tool to discuss the main features of the recon- 
structed f{R)- 

First, let us note that, for i? >> Rq (that is, for very 
small t), the second term in the square brackets may be 
neglected and f{R) ^ R + 771 i?^-^^ Ini?. This model can 
be particularly interesting in the early universe, where it 
allows infiation js^l- In our case, however, this regime 
is achieved also during the recombination epoch, when 
nucleosynthesis took place. In this limit, one should ex- 
pect that standard general relativity (i.e., f{R) oc R) 
is recovered in order to not alter the successful theory 
of primordial nucleosynthesis. Actually, the formation of 
light elements depends on the physics of nuclear reactions 
that, of course, is independent of the form of f{R) and 



Defining e = [f{R) - fapp{R)] /f{R), with / and fapp the nu- 
merically reconstructed and the approximating f{R), we get 
|e| < 0.06 over the full time range 10^''' 1^ t < I, which is 
indeed quite satisfactory. 
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on the expansion rate H{t) setting the beginning and the 
end of the nucleosynthesis process. This latter ingredient 
is the same in our model as in the standard recipe, since 
f{R) has been reconstructed in such a way that H{t) 
has the same functional dependence on t as in the model 
with a scalar field and no departures from general rela- 
tivity. As such, we conclude that having a nonstandard 
f{R) during the recombination epoch does not affect the 
process of formation of primordial elements. 

As the time increases, the f{R) never approaches a 
linear function, so that the general relativity Lagrangian 
f{R) ~ i? is never approached. In particular, for 0.05 < 
T < 0.5 (i.e., I < z < 8), both terms in Eq.(33) con- 
tribute to the Lagrangian giving rise to large departures 
from the standard model. Considering that the typical 
galaxy formation redshift is of order Zi? ~ 2 5, it could 
thus be interesting to investigate whether this could have 
any impact on the theory of galaxy formation. 

For 0.5 < T < 1 (i.e., for z < 0.8, part of the redshift 
range probed by the SNela Hubble diagram), the best 
fitting function significantly differs from standard general 
relativity as expected. Moreover, we note that f{R) is 
neither a simple power law in R (or a combination of 
power law functions), nor a logarithmic one. This could 
be surprising since both these classes of Lagrangians have 
been widely discussed in literature and found to be in 
agreement with the data. Actually, it is well known that 
the data do not select a unique class of models (dark 
energy -like or higher order gravity theories). Moreover, 
the reconstructed f{R) predicts the same dynamics of the 
scalar field model with the exponential potential, so that 
we are sure that it fits the data equally well. However, it 
is intriguing to note that the approximating function is a 
combination of both (power law and logarithmic) classes 
of ,f{R) theories. 



VI. CONCLUSIONS 

Summarizing, we may conclude that the f{R) theory 
reconstructed starting from the exponential potential is 
able to solve the puzzle of cosmic acceleration resorting 
to corrections to general relativity that are significative 
both in the early and late universe. As we have discussed, 
however, considering f{R) oc R during the nucleosynthe- 
sis epoch is not mandatory to get the correct abundance 
of light nuclei, so that this argument may not be used to 
reject the reconstructed /(i?)'s. 

On the other hand, the impact of strong departures 
from general relativity during the galaxy formation epoch 
should be investigated with much care. Actually, while 
the time scale of the process is unaltered since the ex- 
pansion rate of the universe is the same, the physics of 
the gravitational collapse giving rise to the protogalactic 
clouds may be different in a fourth order theory with re- 
spect to quintessence models (independent of the shape 
of the potential). Indeed, the Newtonian potential ob- 
tained considering the weak field limit of a particular 



class of these theories (namely f{R) oc i?") is modified 
and a corrective term as r" has to be added to the 1/r 
Newtonian term [23|- However, it is worth noting that 
it is indeed possible that the strength of the correction 
is too low to significantly alter the process, so that any 
difference could be hardly detected. 

There is a subtle point to be discussed inherently to 
the the peculiar nature of the reconstruction procedure 
we have employed. As explained in Sect. HI, we have set 
boundary conditions at the present time and integrated 
Ea. (|13|l backward in time. As an alternative procedure, 
one could give boundary conditions at some very far point 
in the past and then integrate forward. Indeed, a pos- 
sible strategy could be to impose that f{Ri„) ^ R (and 
hence /'(i?^„) = 1, f"{R^n) = 0) with i?„ = i?(r„) 
and Tin ~ I0~^, that is, to impose that standard general 
relativity is recovered during the recombination epoch. 
Performing such a reconstruction, we infact get a f{R) 
function that departs from the general relativistic one 
{f{R) ^ R) only in the recent epochs, so that all the suc- 
cessful results of the standard theory are automatically 
retained. Although such a feature is highly desirable, this 
approach is manifestly biased by theoretical prejudices. 
After all, there is no aprioristic reason why f{R) should 
be linear in R on scales that have never been directly 
probed by observations. 

It is also interesting to stress that the reconstructed 
f{R) is by no way unique since different boundary condi- 
tions lead to diverse f{R) theories, all sharing the feature 
of giving rise to the same expressions for the Hubble pa- 
rameter and the scale factor. This result does not come 
as a surprise. For each f{R), there is a different La- 
grangian giving rise to a particular set of field equations. 
Since these latter are nonlinear, there may be, in gen- 
eral, more than one solution, and it is thus not unlikely 
that a given solution is common to more than one f{R) 
theory. Going backward from solutions to f{R) there- 
fore cannot uniquely determine f{R), and this explains 
why using different boundary conditions leads to differ- 
ent higher order theories. We stress, however, that this 
cannot be considered as a weakness of our reconstruction 
procedure. Our aim is, in fact, to find out a class of f{R) 
theories able to reproduce a given cosmological dynam- 
ics, and this is indeed what we do successfully. Choosing 
among the possible theories can only be done by care- 
fully investigating what are the most suitable boundary 
conditions. This is outside the aim of the present work, 
but will be discussed in a forthcoming paper. 

While the reconstructed f{R) correctly works on cos- 
mological scales, its consequences on the Solar system dy- 
namics are worth investigating. Actually, this is an open 
problem for all the class of fourth order theories and, 
up to now, no definitive conclusions have been drawn. 
A possible way out is to evaluate PPN parameters for 
the case of f{R) gravity exploiting the analogy between 
fourth order theories and scalar - tensor ones. Such an 
analysis has been indeed performed by some of us [35| 
and will be presented in a forthcoming work. 
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As a concluding remark, we would like to emphasize 
the resulting substantial equivalence between two radi- 
cally different scenarios, namely scalar field quintessence 
and fourth order theory of gravity. Since it is possible 
to find f{R) in such a way that the same dynamics is 
predicted by the two models, infact, we are in presence 
of an intrinsic degeneracy that cannot be broken by ob- 
servations probing only dynamical quantities (but see the 
discussion in the concluding section of j^J). While this is 
bad news for observational cosmology, it is a call to army 
for theoreticians to look for suitable and self consistent 



motivations to choose between these two equivalent ap- 
proaches to the dark energy puzzle. 
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